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We propose a method for the lattice QCD computation of nucleon-nucleon low-energy interactions. 
It consists in simulating QCD in the background of a "electromagnetic" field whose potential is non- 
vanishing, but whose field strength is zero. By tuning the background field, phase-shifts at any (but 
small) momenta can be determined by measuring the shift of the ground state energy. Lattice sizes 
as small as 5 Fermi can be sufficient for the calculation of phase shifts up to momenta of order of 
m-n/2. 



One of the central goals of nuclear physics is to relate the successful phenomenological models developed throughout 
the years with the underlying fundamental theory of the strong interactions, QCD. Effective field theories are an 
important step in this direction, but they are inherently limited by the existence of low energy constants whose 
values are not determined by symmetries and have to be fit to experiment. The need is then obvious for a fully non- 
perturbative method that can determine the interaction between nucleons (or alternatively, the low energy constants 
of the effective theory) directly from QCD. At present, lattice QCD is the only such method. 

Most phenomenological models of nuclei arc based on non-rclativistic two(and three) nuclcon potentials. However, 
since nucleons are not infinitely heavy, the inter-nucleon potential is not a well defined quantity that can be measured 
on the lattice, even in principle. Instead, the connection between QCD and nuclear physics should be established 
through observables like scattering amplitudes and phase shifts, etc.. That brings out a problem: lattice calculations 
are done in euclidean space and analytic continuation of the euclidean correlation functions at infinite volume to 
Minkowski space is, in practice, impossible. This observation, formalized in 0], seems to restrict lattice QCD to 
observables like masses, decays constants and amplitudes at kinematical thresholds. Phase shifts at some special 
values of the momenta can however be obtained by measuring the shifts in the low lying two-particle states due to the 
finite volume H, Q , as long as the lattice size L is larger than the pion Compton wavelength (up to corrections of 
order e"™-^). This can be intuitively understood by realizing that the baryon number two sector of QCD at momenta 
smaller than the pion mass reduces to a non relativistic quantum mechanical system with two nucleons interacting 
through contact interactions. At momenta Q much smaller than the ~ 1/a, where a is the nucleon-nucleon scattering 
length, this contact interaction is perturbative but it becomes strong at Q ~ 1/a. In particular, for lattices with size 
L much larger than the scattering length a the low lying states have typical momenta Q satisfying Q << 1/a, and 
Luscher derived the formula relating the shifts in the energy levels and a as an expansion in powers of a/L. This 
method has been used to obtain pion-pion scattering phase shifts la| but in the two nuclcon sector I am aware of only 
only one quenched calculation performed with a large pion mass Q . 

In the two-nuclcon case the condition L ^ a can hardly be satisfied since the scattering lengths between two 
nucleons are large by QCD standards (5.42 fm in the spin triplet and 23.7 fm in the spin singlet channel) and 
numerical simulations with lattice sizes much larger than this are impractical. For L = a the shifts in the energy 
levels due to the nucleon-nucleon interactions arc not small but can still be reliably computed and used to obtain 
information on the nucleon-nucleon interactions. An analysis of this method in the two-nucleon case was presented 
in 01 . There we found that, after taking into account the strong nucleon-nucleon interactions, lattice sizes L ~ 8 fm 
are necessary for the ground state to have small enough energy for the method to be valid, and even larger sizes if 
the excited states are considered. From the shift in the ground state energy the phase shifts at only one kinematical 
point can be determined. More handles on the phase shifts coming from the excited states would require even larger 
lattice sizes. 

This paper proposes a method that i) allows for smaller lattice sizes and ii) provides information about the phase 
shifts at any momenta smaller than m,r/2. The basic idea is very simple: one just simulates the baryon number two 
sector of QCD in a finite torus and in the background of a fictitious "magnetic" potential with zero field strength, the 
kind of field generated by a thin solenoid going around inside the torus. Due to the Aharonov-Bohm effect |3, lajthe 
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FIG. 1: The lattice with periodic boundary conditions (and two dimensions suppressed) is represented by the surface of the 
outer torus. The fictitious solenoid (inner ring) generates a magnetic vector potential A along direction z (wrapped around the 
torus). The magnetic field is confined inside the ring and vanishes at the surface of the torus, where the lattice is. 



energy levels are changed by this potential despite the fact that the field strength vanishes everywhere on the lattice. 
The strength of the potential can then be adjusted in order to have the ground state to have any energy desired. 
Alternatively we can describe the method as simulating QCD with twisted boundary conditions for the quark in one 
chosen spatial direction. The two descriptions are related by a change of variables amounting to a discontinuous gauge 
transformation. 



A. QCD with the background field 



We will consider QCD in the presence of a C/(l) background gauge field coupling to baryon number of the form 



where z is the unit vector in the z direction and </> is real. 

In the case of two degenerate flavors of Wilson fermions the quark action is 



m„b + 4r x ^ _ 

1 Z^^-x^x 

x./i 



(2) 
(3) 



where b is the lattice spacing,? indexes the lattice sites, fi the directions of the links, r is the Wilson term coefficient, 
Mij defines the quark operator and a sum over flavors is implicit. The link operators f2p(x) are a product of a SU{3) 
matrix C/^(x) and a phase determined by the background field i^fj,{x) = Uf^{yi)e^''^^'^^\ Even in the presence o the 
background field the determinant of M is positive. For that notice that the matrix AI satisfies jc^M^y'y^ = My^ 
, where the dagger means hermitian conjugation on the spin and colors indices only. This relation implies that 
det{M)* = det{M'^) = det^j^Mj^) = det{M). The quark determinant, being a product of a up quark determinant 
and a down quark determinant, is then positive even in the presence of the background field and standard Monte 
Carlo techniques are available. 

Instead of using the quark fields above, satisfying periodic boundary conditions, we can use instead 



= e 
Qx = e' 



■'^ gx 
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(4) 
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FIG. 2: Sum of graphs determining the twO-nucleon scattering amplitude in the effective theory. The vertices include interactions 
with an arbitrary number of derivatives. 

The (7, q fields satisfy twisted boundary conditions at z ^ L: 

qz=N = e~'iqz=o, (5) 

where N is the number of sites in the z direction. 

We will extract the nucleon-nucleon phase shifts from the long (euclidean) time behavior of the finite volume 
correlator 



C{t,p) = (0|T7V^(p,t)P7V(-p,t) . iVt(-p,0)PtiV*(p,0)|0) 

-^e-^°*|(i^;o|7V^(-p,0)ptiV*(p,0)|0)|^ (6) 

t — 'OO 

where P is the projector on the spin triplet or spin singlet channels, N(\i, t) are operators with the quantum numbers 
of nucleons with momentum k at time t and Eq is the energy of the ground state with the quantum numbers of the 
N^-p,t)P^N*{k,t)\0) state. 



B. The effective theory 



For small momenta Q < uit^, the nucleon-nucleon interaction can be described by an effective field theory containing 
only nucleons as explicit degrees of freedom. This effective theory has been used extensively in the computation of 
few-nucleon observables and has been reviewed in |lOl IllL IT^ . It contains only contact interactions, with increasing 
number of derivatives, and we will denote the coefficient of terms with 2n derivatives by C2n- In the presence of the 
background field all derivatives in the effective lagrangian are substituted by covariant derivatives. That is the only 
way that the background field can enter in the effective theory. Terms that are gauge invariant by themselves, for 
instance, anomalous magnetic terms, vanish since the magnetic field vanishes. We can perform a similar change of 
variables as above and eliminate the background field by working with fields satisfying twisted boundary conditions 

V'(x,y,L) = e'*7A(x,y,0) 

i;\x,y,L) = e-'V(a;,y,0) (7) 

For momenta Q <^ 1/a the effective theory is perturbative but for Q ~ 1/a it is non-perturbative. In fact, the 
two-nucleon scattering amplitude is given by the infinite series of diagrams shown in Fig. Let us first consider 
the case of the spin sin glet channel. Using dimensional regularization the diagrams in Fig. ^] can be computed and 
summed. The result is jlolElIT^ 



~ M kcotS-ik ^ ' 

where k is the center-of-mass momentum of the colliding particles, ^ is the renormalization point, S the phase shift 
at momenta fcand the loop sum Iq is^ 



^ Besides the minimal sub traction scheme used here, other schemes were proposed that make the estimates of C2n and power counting 
much simpler lTol llll Il3ll . 
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= -^4^V^' (9) 

where D is the number of space dimensions and we used the standard relation between the amplitude and the phase 
shifts in the second line of Equation (jS)). We can now go to the on-shell point ME = fc^ and relate the phase shifts 
to the constants C2n 

An 1 

kcOtS^ —. (10) 

The same combination of low energy constants appearing above also determine the position of the energy eigenstates 
on a finite volume. They are given by the poles of the finite volume, real time correlator analogue to Eq.© 

C{E,p) = J dte''^^0\TN'^{p,t)PN{-p,t) . N^-p,0)P^ N* {p,0)\0) 

-!(£;„ |ivt(-p,o)ptiv*(p,o)|o)p, (11) 



E 



E-En+iO' 



where, in the spin singlet case, ~ a2T^T^ (the tr's (r's act on spin(isospin)). The computation of this correlator 
in the effective theory receives two kinds of contributions: from s-wave interactions and from higher partial wave 
operators. The first kind are the only ones that survive in the infinite volume limit, if one is careful to either use 
p = or to average the sink and the source in Equation (|ll|l ovcr all possible directions of p. In a finite volume 
and with the background field, p cannot be zero and cannot be averaged over all directions so higher partial wave 
interactions, starting with the p-wave, contribute to Eauation Hll|l . These contribution are suppressed by a factor 
{Q/m-^y^ compared to the leading interactions. Furthermore, for = 0,7r this contributions are further suppressed. 
For (/) = the cubic symmetry forbids the contamination from p-waves. For (j) = it there is an extra "parity" 
symmetry along the z axis that, combined with the two dimensional cubic group in the x — y plane also forbids p-wave 
contributions. As we will see below, values of the background field around « tt are the most interesting ones, so we 
will disregard the higher partial wave pieces in the following. 

The correlator in Equation Hll|l can be computed in the effective theory and the result is (up to corrections of order 
(g/m^)3(0-^)2)) 

where the sum is over all the allowed momenta in the box. For the twisted boundary conditions these allowed momenta 
are 

with nxyTiy and integers. The poles of Equation (|12|l are then determined by 



kLcotS = -Si^^,(b), (14) 
TT 47r^ 



with fc^ = ME and the function S{ri,(j)) is defined by 

S(r],(j))= lim y 47riV (15) 

|ff|<Ar— ' 
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and n = {ux, ny, Uz + 4)/2tt). Notice that S{7], tt + (f>) = S{ri, tt — (f>). In FigurcfBl wc show S'(ry, (j>) as a function of i] 
for a few values of (/>. ^ 

For large values of L, kL cot 6 is also large and Equation p(J|l will be satisfied for values of k where S is close 
to one of its poles. These poles are located at = 47r2(n^ + riy + (n^ + (j)/2Tr)'^)/ L'^ = ^^/i^, {At:'^ + 4)'^)/L'^, • • • , 
corresponding to the eigenstates of free particles in the presence of the background field. Close to the first pole, for 
instance, the function S is dominated by the n = term and we have 

kL cot d ~ 

where (X){4>) is 



kL^otS^- ^^ .... +'^o(0) + O(( ) ), (16) 



co(0) = lim ^2 A ^ - 4^^^ (17) 

n|<JV 

fi#0 



The equation above can be solved iterativcly determining the energy level : 



E 



(j)^ 47r 

Jn? ^ kMi? 



tan 5 



kL 



tan 5 



(18) 



For k <C TO^, kcotS is well approximated by kcotS ^ — 1/a + rQk'^/2 + • • • (effective range expansion), where a is 
the scattering length and tq the effective range. For L ^ ^/aro, fccot^ ^ — 1/a and the formula above reduces to 



E = 



ML^ ^ ML^ 



1 - 



co(0)a 



(19) 



which is the analogue of the "Luscher's formula" 0, 0, 0| • 

For smaller values of L simple expansions as Equation H19|l are not available. Still, we can numeri cally compute 
the function 5(77, (j)) and related energy level in the box with the values of phase shifts. In Figure |4(a)| we show 
the estimate of the ground state energy for boxes of different sizes and for different values of the background field. 
For these estimates we took the values of the phase shifts as given by th effective range formula with parameters 
as = 23.7, ros = 2.73 for the singlet. The validity of our approach is limited to momenta smaller than about half the 
pion mass, corresponding to an energy scale of about 6 McV. Energy states with (the absolute value of the) energy 
larger than this are an artifact of the effective theory and will not exist in a lattice QCD simulation. 



^ A C code computing S(r], tj>) in an efficient way can be downloaded from |http : //www-nsdth ■ Ibl ■ gov/ "bedaque/ 
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(a) singlet (b) triplet 

FIG. 4: Ground state energy (m MeV) of two nucleons in the spin singlet (left) and triplet (right) channel as a function of the 
background field. The three curves correspond, from bottom to top, L — 5,7 and 10 fm. 



Things are a little more complex in the triplet s-wave channel ('^S'l) due to the mixing with the triplet d-wave 
channel {^Di). In the infinite volmxie limit this mixing is generated by the tensor force and is suppressed at low 
energies by a factor of (Q/m^)^ ^M- Two insertions of the tensor force are necessary: one leading from the s-wave 
to the d-wave and another back to the s-wave. Each one of these transitions is of order {Q/m„-)'^ and the total effect 
is ~ {Q/iTiTr)'^. In a finite volume, due to the breaking of rotational symmetry by the shape of the torus and the 
background field, tensor forces can contribute already at leading order. However, it is easy to see by an explicit 
calculation that the contribution linear in the tensor force is proportional to q^q^ — q^S'^^ /i, where q is either an 
internal or external momentum q = p. This contribution vanishes after averaging over spin polarizations. In other 
words, the spin averaged correlator 



C(i.P) = J E {0\TN^{p,t)P'N{-p,t) N^-p,Q)P'^N*ip,0)\Q), (20) 



i=l,2,3 



(P* ~ a a^T2 is the spin triplet projector) will receive contributions to the tensor force only at second order or higher. 
Since the contributions of the tensor force at finite or infinite volume are small (~ (Q/'tItt)^) we will disregard them 
here. 

In Figure 4(b)| we show the estimate of the triplet ground state energy for boxes of different sizes and for different 



values of the background field. We use the values at = 5.425 and rot = 1.75. Again we see that even boxes as small as 
L = 5 f m can support a state with energy small enough to be useful in the extraction of phase shifts, if (j) is adjusted 
to be around (j) ~ tt. In the absence of the background field box sizes of at least L « 8 — 10 fm would be required. 

In both the spin singlet and triplet cases the lower bound on L saturates the minimum value set by finite pion 
mass effects. These effects are suppressed by the factor g-™"^ which for L = 5 fm gives e^'"''^ w 0.03. To use still 
smaller boxes a similar calculation to the one presented here must be done using an effective theory which includes 
pions explicitly and valid for Q ^ <^ rrip. One inconvenient is that calculations in this effective theory will be 
neccesarily truncated to a certain order in the low energy expansion, unlike here where an all orders computation was 
possible. That means that the matching with the lattice results will be more properly described as a computation 
of the effective theory low energy constants than as a determination of the phase shifts. By the other hand this 
calculation can be used to extrapolate to realistic values of the quark mass [la, ^3, UB, ■ 

Besides the possibility of using smaller lattice sizes, another advantage of the background field method described 
here is that the phase shifts at arbitrary kinematical points are probed (as long as they are in the regime k < to^/2). 

The method described here may also be applied to other hadronic interactions. For instance, pion-pion phase 
shifts can be calculated at arbitrary kinematical points by adding a background field coupling to some fiavor charge, 
for instance, electric charge. This spoils the positivity of the determinant so seems to be feasible only in quenched 
calculations. Similarly, it might be of use in ameliorating some of the issues involved in the k — > tttt lattice extraction. 
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